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Abstract 

The three body decays B — >■ D^*^ D^*'^ Ks may be used to measure both sin 2/? and 
cos 2/3. Crucial to the cos 2/3 measurement is the resonant contribution to the three 
body decay from p-wave excited Ds states. If these p-wave states are the newly 
discovered Ds{2317) and Ds(2460) then they are below the D^*^K threshold and 
hence do not contribute to 5 — > D^*'^ D^*^ Ks- The three body decays can then be 
used to measure sin 2/3 without resonant dilution and to look for new physics in 
b —>■ CCS transition. 
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1 Introduction 



The decay J/ipK^ provides a clean measurement of the angle sin(2/?) in the 

unitarity triangle^]. Both BaBar and Belle have measured this CP phase, with the 
world average being j2] 

sin 2/? = 0.736 ±0.049 . (1) 

Other modes can also provide relevant information on the angle (3, an example 
being the decay — >• D'^*'^D^*\ The possibility of extracting cos 2/? from the decay 
B^ DDKs was mentioned in Ref. P]. These modes are enhanced relative to 

5° D^*W^*^ by the factor ~ 20. As in the case of S° ^ J/tpK^ decay, 

the penguin contamination is expected to be small in these decays. Moreover these 
decays can be used to probe both sin 2/5 and cos 2/5 which can resolve /3 — > 7r/2 — /3 
ambiguity flj. 

2 (] from B D^*W^*^Ks 

The amplitude for the decay 5° D*D*Ks can have a resonant contribution and a 
non-resonant contribution. For the resonant contribution the D*Ks in the final state 
comes dominantly from an excited Ds{l'^) state. In the approximation of treating 
D*D*Ks as D* D s{excited) , there are four possible excited p-wave Dg states which 
might contribute. These are the two states with the light degrees of freedom in 
a = 3/2+ state and the two states with light degrees of freedom in a = 
1/2+ state. Since the states with = 3/2+ decay via d-wave to D*Ks, they are 
suppressed. Of the states with light degrees of freedom in j"^ = 1/2+ states, only 
the 1+ state contributes. The 0+ state is forbidden to decay to the final state D*Ks. 

To estimate the above contribution and to calculate the non-resonant amplitude, 
we use heavy hadron chiral perturbation theory (HHCHPT)|n]. The momentum 
of Kg can have a maximum value of about 1 GeV for B^ D*^ D*~ Kg. This is 
of the same order as which sets the scale below which we expect HHCHPT to 
be valid. It follows that in the present case it is reasonable to apply HHCHPT to 
calculate the three body decays. 

In the lowest order in the HHCHPT expansion, contributions to the decay am- 
plitude come from the contact interaction terms and the pole diagrams which give 
rise to the non-resonant and resonant contributions respectively. The pole diagrams 
get contributions from the various multiplets involving Dg type resonances as men- 
tioned above. In the framework of HHCHPT, the ground state heavy meson has 
the light degrees of freedom in a spin-parity state = \ , corresponding to the 
usual pseudoscalar- vector meson doublet with = (0~, 1~). The first excited state 
involves a p-wave excitation, in which the light degrees of freedom have = | or 
1^. In the latter case we have a heavy doublet with = (l+,2+). These states 
can probably be identified with Dgi{2536) and Dgj{2573). Heavy quark symmetry 
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rules out any pseudoscalar coupling of this doublet to the ground state at lowest 
order in the chiral expansion |H] ; hence the effects of these states will be suppressed 
and we will ignore them in our analysis. 

The other excited doublet has = (0+, 1+). These states are expected to decay 
rapidly through s-wave pion emission and have large widths [7J. The l"*" state in the 
D system has already been seen. Only the 1"^ can contribute in this case. For later 
reference, we denote this state by D*^. However, quark model estimates suggest jH] 
that these states should have masses near m + 6m with 6m = 500 MeV, where m is 
the mass of the lowest multiplet. 

Assuming that the leading order terms in HHCHPT give the dominant con- 
tribution to the decay amplitude we will neglect all sub-leading effects suppressed 
by and 1/m, where m is the heavy quark mass. It can be shown that from 

the time dependent analysis of D*^ D*^ Kg one can extract sin(2/5) and 

cos(2/5). Measurement of both sin(2/5) and cos(2/3) can resolve the [3 7i/2 — (3 
ambiguity as already mentioned. The measurement of sin(2/?) can be made from 
the time dependent partial rate asymmetry while a fit to the time dependent rate 
for T[B^{t) D*+D*'Ks] + T[B%t) D*+D*~K,] may be used for the extraction 
of cos(2/5). The cos(2/9) term measures the overlap of the imaginary part of the 
amplitudes for B D*~^D*~Ks and B — * D*^ D*^ Kg decays and is non zero only 
if there is a resonance contribution. 

As in the case for B D*~^D*~ the asymmetry in B ^ D*^D*^Ks is also 
diluted. For the non resonant contribution to B ^ D*^D*~Ks the final state is 
an admixture of CP states with different CP parities. This leads to the dilution of 
the asymmetry and this is the same dilution of the asymmetry as in the case for 
B — > D*~^D*~ . When the resonant contribution is included there is further dilution 
of the asymmetry from the additional mismatch of the amplitudes for B and B 
decays. One can reduce the additional dilution of the CP asymmetry by imposing 
cuts to remove the resonance. A narrow resonance is preferable as it can be more 
effectively removed from the signal region than a broad resonance . When we include 
the resonance contribution it turns out that a broader resonance leads to a larger 
value of D and is a more useful probe of cos(2/?) because of the the larger overlap 
of the amplitudes for B D*+D*-Ks and B D*+D*-Ks decays. 

The differential decay distribution of the time independent process B^ D*^D*^ Kg 
can be used to discover the l"*" resonance D*^^ if it is above the D^*'iK threshold. The 
differential decay distribution for small values of E^, the kaon energy, shows a clear 
resonant structure which comes from the pole contribution to the amplitude with 
the excited = 1"*" intermediate state. Therefore, examination of the D*Ks mass 
spectrum may be the best experimental way to find the broad 1+ p-wave Dg meson 
and a fit to the decay distribution will measure its mass and the coupling. 

The extraction of sin 2/5 and cos 2/5 from the time dependent rate for B{t) 
D*^D*~Ks can be done in the following manner: We define the following amplitudes 

^MM ^ A{B\p) ^ Dl:{p^)D-^:{p.)Kg{pu)) (2) 
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^A,A. ^ A{B'{p)^Dt:{p^)D-,:{p_)KM, (3) 

where S° and S° represent unmixed neutral B and Ai and A2 are the polarization 

indices of the D*^ and D*~ respectively. 

The time-dependent amplitudes for an oscillating state B^[t) which has been 
tagged as a B^ meson at time t = is given by, 

^A„A.(^) ^ a^i.A.cos(^)+ie-^^V^'^^sin(^) , (4) 

and the time- dependent amplitude squared summed over polarizations and inte- 
grated over the phase space angles is: 

\A{s^,s~;t)\^ = ^ [Go(s+,s") + Gc(s+,s")cosAmit-Gs(s+,s")sinAmt 

Zi 

with 

Go(s+,s") = |a(s+,s-)p+ |a(s+,s-)p, 

Ge(s+,s-) = \a{s+,s-)\'-\a{s+,s-)\', 

G,{s^,s-) = 2^5 (e-2^^a(s+,s-)a*(s+,s-)) 

= -2 sin(2/3) 3? (aa*) + 2 cos(2/3) S (aa*) . 

The variables s+ and s~ are the Dalitz plot variable 

s+ = {p+ + Pfc)^ = (p_ + Pfc)^ 

The transformation defining the CP-conjugate channel B^{t) D*~D*^Ks is s"*" 
s~, a ^ a and /3 — > — /3. Then: 

= ^ [Go(s",s+) -Gc(s",s+)cosAmt-FGs(s",s+)sinAmi 

Zi 

Note that for simplicity the e"^* and constant phase space factors have been omitted 
in the above equations. 

It is convenient in our case to replace the variables and s~ by the variables 
y and where E^ is the Kg energy in the rest frame of the B and y — cos 9 with 
9 being the angle between the momentum of Kg and D*'^ in a frame where the two 
D* are moving back to back along the z- axis. This frame is boosted with respect 
to the rest frame of the B with f3 = —{pk/mB){l/{l — Ek/mB))- Note s+ ^ s~ 
corresponds to y ^ —y. The variable y can be expressed in terms of variables in 
the rest frame of B. For instance 

ms 
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where and E'j^ are the energy of the D*^ in the rest frame of the B and in the 
boosted frame while E'^ is the energy of the B in the boosted frame. The magnitudes 
of the momentum of the B and the D*~^ in the boosted frame are given by p'^ and 
p'_^ respectively. 

In the approximation of neglecting the penguin contributions, proportional to 
the small CKM elements VubV*^, to the amplitude there is no direct CP violation. 
This leads to the relation 

a^''^^{pki,Ek) = a-^''-^'{-p„„E,) (6) 

where pki is the momentum of the of the Kg in the boosted frame. The above 
relations then leads to 



Goi-y,Ek) = Go{y,Ek) (7) 

G,i-y,Ek) = -G,{y,Ek) (8) 

Gsi{-y,Ek) = Gsi{y,Ek) (9) 

Gs2{-y,Ek) = -G,2{y,Ek) (10) 



where we have defined 



Gsi{y,Ek) = ^{aa*) (11) 
Gs2{-y,Ek) = Q{aa*) (12) 

Carrying out the integration over the phase space variables y and E^ one gets the 
following expressions for the time-dependent total rates for B^[t) — *• D*~^D*~Kg and 
the CP conjugate process 

r(t) = ^[/o + 2sin(2/3)sin(Amt)/,i] (13) 
r(t) = ^[/o-2sin(2/3)sin(Amt)/,i] (14) 

where Iq and Igi are the integrated Go{y,Ek) and Gsi{y,Ek) functions. One can 
then extract sin(2/5) from the rate asymmetry 

^[^j ~Ej^j = D sin(2/5) sin(Amt) (15) 

where 

D='^ (16) 
is the dilution factor. The quantities J^i o can be calculated in HHCHPT^. 
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The cos(2/9) term can be probed by by integrating over half the range of the 
variable y which can be taken for instance to be y > 0. In this case we have 

V{t) = i[Jo + ^cCos(Amt) + 2sin(2/5) sin(Amt)J^i - 2cos(2/5) sin(Amt)J,2] 

T{t) = ^[Jo + JcCos{Amt) - 2sin(2/3) sin(Amt)J^i - 2cos(2/3) sin(Amt)J,2] 

where Jo, Jc, Jsi and Js2, are the integrated Go{y,Ek), Gc{y,Ek), Gsi{y,Ek) and 
Gs2{y,Ek) functions integrated over the range y > 0. The quantities Jo,c,si,s2 can 
be calculated in HHCHPTj4j. One can measure cos(2/5) by fitting to the time 
distribution of T(t) + T(t). Measurement of the cos(2/?) then resolves the /? — | — /3 
ambiguity. In passing we note that only the sign of cos(2/3) is required to resolve 
the the /? — * | — /? ambiguity. 

3 The new Dgj resonances 

The previous year saw the discovery of an unexpectedly light narrow resonance in 
D+TT^ with a mass of 2317MeV/c'^ first reported by the BaBar collaboration [Oj, 
together with another second narrow resonance in Dgir^'y with a mass 2460MeV^/c^ 

m- 

The smaller than expected masses and narrow widths of these states have led, 
among other explanations to a multi-quark anti-quark or a DK molecule in- 
terpretation of these states |T2], or to an interpretation as p-wave states where the 
light degrees of freedom are in an angular momentum state jq = | or even 
some combination of these ^Ij. There are also conflicting lattice interpretations of 
these states [15 . The mass difference between the Ds{2317) and the well estab- 
lished lightest charm-strange meson, Dg, is AM = 350MeV/c^. This is less than 
the kaon mass, thus kinematically forbidding the decay Ds(2317) D^^d + K. The 
possible resonance at 2460Me V"/ also has such a mass difference when taken with 
the lighter D* state. The interpretation of these states as bound D^*^K molecules 
just below the D^*^K threshold is particularly interesting in the light of the recent 
discovery of a narrow resonance in the decay J/ip —>■ jpp ^H] which has been inter- 
preted as a zero baryon number, "deuteron-like singlet ^Sq" bound state of p and p 

m- 

The production of these new D^j states in non leptonic B decays appear to be in 
conflict with theory predictions based on factorization and heavy quark symmetry 
[TH| ITU] . Following Ref . [18^ let us first assume that we can identify the the newly 
discovered states £',(2317) with D^o and £',(2460) with D*^. In the Standard Model 
(SM) the amplitudes for B D^*^ Dso(D*-^), are generated by the following effective 
Hamiltonian j20j: 
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E(K6K:,cr + Vay:,c'i + Vt,v:^c'M] + h.c. , (i?) 



i=3 



where the superscript u, c, t indicates the internal quark, / can be u or c quark, 
q can be either a ci or a s quark depending on whether the decay is a AS* = or 
AS' = — 1 process. The operators Of are defined as 



= q^j^Lfpfp^Lb^ , Olf = qj,LffrLb, 

01,5 = qif^Lbq'-ff,L{R)q' , = qa7^Lbf3q'f^-ff,L{R)q'^ , (18) 

3 3 
^7,9 = ^qif.Lbeg,q'YR{L)q' , OI^q = -qa-ii,Lbpeq>q'p-f^R{L)q'^ , 

where R{L) = 1 ± 75, and q' is summed over all fiavors except t. O1/2/ are the 
current-current operators that represent tree level processes. O^-q are the strong 
gluon induced penguin operators, and operators O7-1Q are due to 7 and Z exchange 
(electroweak penguins), and "box" diagrams at loop level. The values of the Wilson 
coefficients can be found in Ref. |20j . 

In the factorization assumption the amplitude for B D^*^Dso{D*{), can now 
be written as 

M = M1 + M2 (19) 



where 



Ml = <D,o(^:i)|^7m(1-7')c|0><dW|c7"(1-7')&|5> 
M2 = <Z^.o(/^:i)|s(l + 7')c|0><Z}W|c(l-7^)6|fi> (20) 



V2^ 



(21) 



where 



We have defined 



with 



X2 = -2[ ^B, + B, + -^Br + ^8 ) (22) 



B^ = - Y: 4Vq (23) 

q=u,c,t 



Vq = V;,V,, (24) 

In the above equations Nc represents the number of colors. To simplify matters 
we neglect the small penguin contributions and so as a first approximation we will 
neglect M2. 
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We can now define the following ratios 

BR\B ^ DDso] 



BR[B 




BR[B - 




BR[B - 


D*Ds] 


BR[B - 


- dd:i] 


BR[B 


- DDt] 


BR[B 


D*D*J 



Rdo 
Rd*o - 
Rdi 

- BR[B ^ D*Dt] ^^^^ 

Let us focus on the ratio -Rdo which within factorization and the heavy quark limit 
can be written as 

Rdo = (26) 

where we have neglected phase space ( and other) effects that are subleading in the 
heavy quark expansion. Similarly we have 

Rdi = (27) 

Now in the heavy quark limit f^^o = /d*^ and fn^ = fn* and so one would 
predict Rdq ~ -Rdi- There have been various estimates of the decay constant fc^o 
in quark models and in QCD sum rule calculations ; these typically find the p-wave 
, jq = \ states to have the similar decay constants as the ground state mesons. We 
therefore expect foso ~ Jd^ giving in addition to the heavy quark predictions 

Rdo ^ ~ 1 (28) 

Experimentally Belle measures [22] 

BR[B DDs{2317)]BR[Ds{2317) D,7r°] = (9.9l|^ ± 3.0) x 10"^ 

The dominant decay of the 1^5(2317) is expected to be through the DgTi mode [SB] 
and so 

BR[D ^ DDs{2317)] ^ 10"^ (29) 
Now using the measured branching ratio 



BR[B+ = (1.3 ± 0.4) X 10"^ 

BR[Bd ^ D-D^] = (8 ± 3) X 10~^ (30) 
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one obtains a combined branching ratio 

BR[B DDs] ^ 10^2 (31) 

This leads to -R^o ~ {oyJd^^ ~ |/z)J which is a factor 10 smaller then theoretical 
expectations. 

One might argue that factorization is not applicable to B ^ [)(*)[)(*) decays. 
However recent analysis in Ref. [23] find that factorization works well for these 
decays. Moreover the quantities in Eq. |231are ratios of nonleptonic decay amplitudes 
and so nonfactorizable effects may cancel. So what one really requires is significantly 
different nonfactorizable corrections between decays with the p-wave states in the 
final state and decays with the ground state mesons in the final state. It is possible 
that the discrepancies between experiments and theory may arise from a combination 
of incorrect model prediction of p-wave state properties and nonfactorizable effects. 

If these new Dgj states are indeed the p-wave cs meson then they cannot con- 
tribute to B ^ D^*W^*^Ks as they are below the D^*^K threshold. Hence to the 
extent that the corrections to predictions of HHCHPT are small the three body de- 
cays B D^*^ D^*^ Ks should be dominated by the non resonant contribution. This 
would then imply that sin 2/? can be cleanly measured in this decay. On the other 
hand if the new Dsj states are not the p-wave cs mesons but something exotic like 
four quark states or molecules then experimentally the real p-wave states should 
show up. 

4 New Physics 

There are many reason to believe that the Standard Model is not a complete theory 
as it leaves several puzzles unresolved, specially in the flavour sector. There are 
several hints of possible deviations from the SM [^Hl 123 and several interesting 
methods have been proposed to measure the parameters of the underlying new 
physics j2Zl . If there is new physics in b —>■ s transition then this should show up in 
decays with the underlying quark transition b —>■ ccs which can then be probed in 
B D^*W^*^Ks. One can probe for this new physics also in other decays with the 
same quark level transition such as -B ^ J/ipKs- However, it is possible that the 
matrix element of the new physics may be suppressed in some decays and enhanced 
in other making it mandatory to study as many different decays as possible. If the 
new physics involves a new effective b sg vertex its effect in 5 — > J/tpK^ may 
be negligible because of OZI suppression. Such a suppression would not apply to 
B D^*W^*^Ks decays. 
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5 Conclusion 



In this talk we have pointed out that the three body decays B D^*^ D^*^Ks may be 
used to measure both sin 2/? and cos 2/5. The cos 2/? measurement requires resonant 
contribution to the three body decay from p-wave excited Dg states. We discussed 
the newly discovered Dgj states and pointed out that their properties are sometimes 
inconsistent with an interpretation as p-wave cs states. If these p-wave states are the 
newly discovered Ds(2317) and -0^(2460), , then they are below the D^*'>K threshold 
and hence do not contribute to B ^ D^*'^ D'^*^ Ks- The three body decays can then 
be used to measure sin 2/5 without resonant dilution and to look for new physics in 
b CCS transition. 
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